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Abstract 

In  the  present  work  the  operation  of  a  SOFC  based  on  a  co-ionic  electrolyte  is  considered.  The  special  case  of  ion  transport  in  a  co-ionic 
electrolyte  is  analyzed.  It  is  shown  that  partial  ion  currents  do  not  correspond  to  ion  transfer  numbers,  a  partial  proton  current  being  higher 
than  a  proton  transfer  number  under  the  working  conditions  of  the  SOFC.  It  is  shown  that  both  fuel  utilization  and  electric  efficiency  are 
higher  when  the  proton  transfer  number  is  higher.  It  is  established  that  the  maximum  achievable  efficiency  of  a  hydrogen-fed  SOFC  based  on 
the  co-ionic  electrolyte  with  t\\  =  0.5  is  0.78  at  900  K  whereas  the  efficiency  is  0.62  when  the  SOFC  works  at  70%  of  its  maximum  power. 
©  2004  Elsevier  B.V.  All  rights  reserved. 
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1.  Introduction 

Solid  oxide  fuel  cells  (SOFCs)  are  considered  as  one  of 
the  most  promising  electric  generators.  At  the  initial  stage,  a 
SOFC  based  on  an  oxygen  ion  electrolyte,  SOFC(02~),  has 
been  developed.  During  the  last  two  decades  however,  a  lot 
of  oxides  possessing  proton  conduction  have  been  discov¬ 
ered  [1].  They  were  called  “proton  conductors”  or  “proton 
electrolytes”.  Some  of  them  have  purely  proton  conduction, 
but  the  majority  have  both  proton  and  oxygen  ion  conduc¬ 
tion.  They  can  be  called  “co-ionic  electrolytes”  [2]  and  char¬ 
acterized  by  a  proton  transfer  number,  tu,  which  is  usually 
measured  in  “one-chamber”  conditions  [3]. 

Theoretical  analysis  showed  that  a  SOFC  based  on  a 
pure  proton  electrolyte,  SOFC(H+),  can  reach  considerably 
higher  efficiency  than  the  SOFC(02-)  [4].  It  is  of  theoreti¬ 
cal  and  practical  interest  to  study  the  influence  of  the  proton 
transfer  number  on  efficiency  and  other  characteristics  of 
the  SOFC  based  on  a  co-ionic  electrolyte,  called  hereafter 
SOFCXO2-,  H+). 

2.  Theoretical  model 

A  co-ionic-electrolyte-based  cell  and  its  equivalent  circuit 
are  presented  in  Fig.  1.  Here  it  is  assumed  that  oxidant  is 
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at  the  left  side  and  fuel  is  at  the  right  side  of  the  cell.  It 
is  assumed  that  Ri  =  R/U,  where  R  is  ohmic  resistance 
of  the  electrolyte,  6  is  corresponding  ion  transfer  number. 
Electromotive  forces  corresponding  to  each  kind  of  charge 
carrier  are  calculated  by  equations 


RT  p"  02 

- In - , 

4  F  p'02 

RT  P"H  2 

— In - 

IF  P'  H2 


(1) 
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where  pu2  at  the  left  side  of  the  cell  and  po2  at  the  right  side 
are  calculated  using  the  equilibrium  constant  of  the  reaction 

H2  +  3O2  =  H20  (3) 


The  electromotive  force  of  the  cell  represents  a  stationary 
one  and  is  equal  to 

E  =  tuEu  +  t0E0  (4) 


Current  density  across  the  co-ionic  electrolyte  at  a  certain 
point  depends  on  the  emf  at  this  point  and  the  terminal 
voltage  U : 


E-U 

j  =  - , 

Pc  ff 


where  peff  is  the  total  resistance  of  1  cm2  of  the  cell.  Par¬ 
tial  ion  current  densities  are  determined  by  corresponding 
electromotive  forces  as  follows  from  the  equations 
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Fig.  1.  An  electrochemical  cell  with  a  co-ionic  electrolyte  (left),  its  equivalent  circuit  (center),  and  a  scheme  of  a  fuel  cell  based  on  a  co-ionic  electrolyte 
(right). 
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Generally,  Eft  ^  Ep,  so  the  partial  ion  currents  do  not  cor¬ 
relate  to  the  transfer  numbers.  Moreover,  the  U  value  can  be 
between  Eft  and  Eo  and  one  of  the  partial  currents  can  be 
negative. 

In  the  case  of  a  single  cell  its  efficiency  can  be  calculated 
by  the  following  equation: 

U 

ri  =  m~rr ,  (8) 

Uib 


where  rj f  is  the  fuel  utilization,  f/t h  is  the  thermoneutral 
voltage  defined  as  follows 
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o 
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where  AH ^  is  the  enthalpy  of  the  fuel  oxidation  reaction. 
The  upper  limit  of  the  efficiency  is  determined  by  the  value 
of  the  cell  voltage  at  which  the  total  current  density  is  pos¬ 
itive  in  any  point  of  the  cell  at  a  minimum  admissible  par¬ 
tial  pressure  of  the  fuel  components  in  the  anode  mixture. 
This  voltage  value  can  be  called  the  “maximum  admissible 
voltage”,  or  simply  “maximum  voltage”  of  the  fuel  cell  at 
such  fixed  operating  parameters  as  temperature,  input  fuel 
gas  composition,  input  air  humidity,  minimum  concentration 
of  fuel  components  in  the  anode  output  gas,  and  oxygen  uti¬ 
lization  in  the  cathode  gas  (hereafter  operating  parameters). 
The  maximum  voltage  can  be  calculated  by  the  following 
equation 


a 


max 


fif(l  ~  PH2mm)Eft  +  tpEp 
1  —  fifPH2min 


(10) 


where  PH2min  is  the  minimum  admissible  hydrogen  partial 
pressure  in  the  anode  gas.  On  the  whole,  the  maximum  ter¬ 
minal  voltage  is  slightly  lower  than  the  cell  emf  at  the  inlet. 

An  important  operating  performance  of  the  fuel  cell  is 
its  relative  power  pr ,  equal  to  the  ratio  of  its  current  fuel 
cell  power,  P,  to  the  maximum  achievable  power  (at  fixed 
operating  parameters)  PmSLX'- 

Pr  =  ~~  ■>  (ID 

max 


The  fuel  cell  efficiency  is  related  to  the  relative  power: 

0  ~  0.5/7max(l  +  yf  1  px) ,  (12) 


where  77max  is  the  upper  achievable  efficiency  at  fixed  op¬ 
erating  parameters.  The  fuel  cell  operation  at  high  pv  when 


its  voltage  is  close  to  half  of  the  emf  is  unacceptable  be¬ 
cause  of  its  low  efficiency  (ca.  50%  of  the  maximum)  and, 
hence,  high  current  losses.  The  fuel  cell  operation  at  low  pY 
is  attractive  from  the  point  of  view  of  high  efficiency  (and, 
hence,  low  current  losses)  but  this  is  reasonable  only  when 
the  cost  of  the  fuel  cell  is  low.  Development  of  SOFC  tech¬ 
nology  leads  to  a  reduction  of  prices  of  the  SOFC  units,  so 
analysis  of  the  SOFC  operation  at  high  efficiency  is  impor¬ 
tant  for  understanding  the  potential  of  SOFCs.  One  of  the 
acceptable  regimes  of  fuel  cell  operation  is  at  pv  ~  0.7.  In 
this  case  the  fuel  cell  power  can  be  easily  changed  (in  par¬ 
ticular,  the  power  can  be  rapidly  increased  by  ca.  40%)  at 
the  same  operating  parameters  and  the  fuel  cell  efficiency 
in  this  case  is  relatively  high  (ca.  77%  of  the  maximum 
efficiency). 


3.  Results  and  discussion 

The  following  operating  parameters  are  accepted  in  the 
frame  of  this  paper:  the  anode  input,  humid  (2%)  hydrogen; 
Pft 2  =  0.1  in  the  anode  output;  the  cathode  input,  humid 
(2%)  air,  the  oxygen  utilization  0.1;  the  fuel  cell  tempera¬ 
ture  900  K.  It  is  assumed  that  tft  =  tp  =  0.5  and  do  not  de¬ 
pend  on  anode  and  cathode  gas  composition.  The  results  are 
presented  for  two  regimes  of  the  SOFC  operation:  at  maxi¬ 
mum  terminal  voltage  and  at  pv  =  0.7,  and  for  two  modes 
of  the  SOFC  feeding:  co-flow  and  counter-flow  modes.  Dis¬ 
tributions  of  parameters  are  presented  for  the  dimensionless 
channel,  i.e.  in  coordinates  v/L,  L  the  channel  length. 

In  order  to  demonstrate  the  special  case  of  the  SOFC(02~ , 
H+)  more  clear  by  the  most  important  characteristics  of  the 
SOFC(02-)  and  of  the  SOFC(H+)  for  the  case  of  the  max¬ 
imum  voltage  are  shown  in  Figs.  2-6.  Charge  transfer  and 
hence  the  electrical  characteristics  of  the  SOFC  are  is  deter¬ 
mined  by  the  distribution  of  reactant  components  along  the 
SOFC  channels.  Gas  mixtures  in  the  fuel  cell  channels  can  be 
described  by  flow-rates  of  the  gas  components  and  by  their 
partial  pressures.  Distributions  of  the  anode  components  par¬ 
tial  pressures  and  their  flow-rates  expressed  as  A^/Ah2(0), 
where  A;  is  flow-rate  of  /-component  and  Ah2  (0)  is  hydro¬ 
gen  flow-rate  at  the  inlet,  for  the  case  of  maximum  voltage 
are  plotted  in  Fig.  2  for  the  SOFC(02-).  It  is  evident  that 
total  gas  flow-rate  is  constant  and  only  redistribution  of  the 
partial  flow-rates  and  of  the  partial  pressures  of  the  compo¬ 
nents  occurs  in  the  anode  channel.  It  is  possible  to  show  that 
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Fig.  2.  Distributions  of  the  anode  components’  relative  flow  rates  ( — ) 
and  the  anode  components’  partial  pressures  ( — )  along  the  S0FC(02-). 
Case  of  maximum  voltage. 


Fig.  5.  Distributions  of  the  emf  along  the  SOFC(H+)  and  S0FC(02-)  for 
the  co-  (thin  lines)  and  counter-flow  (bold  line)  modes.  Case  of  maximum 
voltage. 


Fig.  3.  Distributions  of  the  cathode  components’  partial  pressures  along  the 
S0FC(02-)  for  the  co-flow  mode  (dashed  lines),  and  along  the  SOFC(H+) 
for  the  co-flow  (thin  lines)  and  counter-flow  (bold  lines)  modes. 


the  fuel  utilization  can  be  expressed  as 
Nu2(L) 


'?/  =  !- 


VH,  (0)  ’ 


(13) 


One  can  see  that  the  fuel  utilization  in  the  SOFC(02  )  is 
about  0.9. 
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Fig.  6.  Distributions  of  the  current  densities  along  the  SOFC(H+)  and 
S0FC(02-)  for  the  co-  (thin  lines)  and  counter-flow  (bold  line)  modes. 
Case  of  maximum  voltage. 

Distributions  of  the  cathode  components’  partial  pressures 
are  shown  in  the  Fig.  3.  In  the  S0FC(02-),  the  oxygen 
partial  pressure  changes  slightly.  The  consequence  of  this 
is  an  absence  of  noticeable  dependence  of  the  cathode  po¬ 
tential  on  both  the  channel  coordinate  and  the  mode  of  the 
S0FC(02-)  feeding.  In  the  SOFC(H+),  due  to  significant 
relative  change  of  the  steam  partial  pressure,  the  cathode 
potential  changes  considerably  from  the  cathode  inlet  to  the 
outlet.  This  causes  a  noticeable  difference  in  the  SOFC(H+) 
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Fig.  4.  Distributions  of  the  anode  components’  relative  flow  rates  (solid  lines)  and  the  anode  components’  partial  pressures  (dashed  lines)  along  the 
SOFC(H+)  for  the  cases  of  co-flow  (left)  and  counter-flow  (right)  modes.  Case  of  maximum  voltage. 
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Fig.  7.  Distributions  as  in  the  Fig.  2  along  the  SOFC(02 


,  H+)  for  the  co-  (left)  and  counter-flow  (right)  modes.  Case  of  maximum  voltage. 


characteristics  for  different  feeding  modes.  Fig.  4  illustrates 
the  influence  of  the  feeding  mode  on  the  characteristics  of 
the  anode  flow  in  the  SOFC(H+).  The  figures  also  show 
that  the  total  anode  gas  flow-rate  decreases  from  the  inlet  to 
the  outlet  and  becomes  practically  equal  to  the  steam  par¬ 
tial  flow-rate.  One  can  see  that  the  fuel  utilization  in  the 
SOFC(H+)  is  very  close  to  1  in  the  maximum  voltage  case. 

Fig.  5  demonstrates  distributions  of  the  electromotive 
forces  along  the  two  types  of  the  SOFCs.  As  is  mentioned 
above,  the  emf  distribution  in  the  S0FC(02-)  does  not 
depend  on  the  feeding  mode.  On  the  contrary,  the  emf  dis¬ 
tribution  in  the  SOFC(H+)  depends  strongly  on  this  factor. 
One  of  the  consequences  of  this  is  a  difference  of  the  emf  at 
the  SOFC(H+)  outlet  for  different  feeding  modes:  E(L)  = 
1013  mV  for  the  co-flow  mode  and  E(L )  =  1057  mV  for 
the  counter-flow  mode.  Note  that  E(L)  =  909  mV  for  the 
S0FC(02-).  It  is  clear  that  the  terminal  cell  voltage  cannot 
be  higher  than  E(L ),  so  the  SOFC(H+)  has  an  evident  ad¬ 
vantage  against  the  S0FC(02-)  as  regards  the  possibility  to 
reach  higher  efficiency,  the  counter  flow  mode  being  more 
attractive  for  this  purpose. 

Distributions  of  current  densities  presented  in  the  form 
jyOeff  are  shown  in  Fig.  6.  A  strong  irregularity  of  the  current 
density  in  the  S0FC(02-)  causes  significant  thermal  stresses 
due  to  the  irregularity  of  heat-evolution.  The  current  density 


Fig.  8.  Distributions  of  the  electromotive  forces  corresponding  to  various 
charge  carriers  along  the  SOFC(02-,  H+)  for  the  co-  (thin  lines)  and 
counter-flow  (bold  lines)  modes.  Case  of  maximum  voltage. 


within  most  of  the  SOFC(H+)  changes  considerably  less 
than  in  the  previous  case.  Numerical  integration  allows  us 
to  find  an  average  current  per  unit  of  the  SOFC  area,  yav, 
and  consequently  an  average  specific  SOFC  power,  Pav, 


(£av  -  U)U 
Peff 


(14) 


Figs.  7-9  illustrate  characteristics  of  the  S0FC(02  ,  H+) 
for  the  maximum  voltage  case.  The  total  anode  flow-rate  as 


Fig.  9.  Distributions  of  the  current  densities  along  the  SOFC(Q2  ,  H+)  for  the  co-flow  (left)  and  counter-flow  (right)  modes.  Case  of  maximum  voltage. 
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Fig.  10.  Distributions  as  in  Fig.  2  along  the  SOFC(02  ,  H+).  Case  of 
Pr  =  0.7. 

well  as  the  hydrogen  flow-rate  slow  down  gradually  prac¬ 
tically  to  zero  at  the  end  of  the  channel  while  the  steam 
flow-rate  reaction  reaches  a  maximum  (Fig.  7).  The  hydro¬ 
gen  partial  pressure  reaches  a  minimum  value  at  the  end  of 
the  channel  in  the  case  of  the  counter-flow  mode,  while,  in 
the  case  of  the  co-flow  mode,  it  has  a  minimum  approxi¬ 
mately  in  the  middle  of  the  channel  and  then  rises  slightly. 
It  is  necessary  to  point  out  that  neither  the  operating  regime, 
nor  the  terminal  voltage  affect  noticeably  the  nature  of  the 
distributions  of  the  cathode  components  and  they  are  prac¬ 
tically  the  same  as  is  reported  in  Fig.  3. 

Fig.  8  demonstrates  distributions  of  the  electromotive 
forces  corresponding  to  various  charge  carriers.  The  values 


Fig.  11.  Distributions  as  in  Fig.  8.  Case  of  pY  =  0.7. 

of  the  emf  at  the  end  of  the  cell  differ  considerably.  As  is 
mentioned  above,  the  terminal  voltage  can  be  between  these 
values:  Eq(L)  <  U  <  En(L).  In  this  case,  the  oxygen  ion 
current  density  is  negative  while  the  total  current  density  is 
positive  (see  Fig.  9).  Of  course,  a  substantial  part  of  the  cell 
generates  an  insignificant  part  of  the  cell  power  and  could 
be  eliminated,  but  this  would  lead  to  a  lowering  of  the  fuel 
utilization  and  hence  to  a  lowering  of  the  efficiency.  On  the 
whole,  the  current  density  is  very  irregular,  although  the 
counter-flow  mode  is  more  favourable. 

Characteristics  of  the  S0FC(02-,  H+)  at  pr  =  0.7  are 
shown  in  Figs.  10-12.  An  increase  of  the  relative  power,  i.e. 
a  decrease  of  the  terminal  voltage,  influences  strongly  all  the 
cell  characteristics.  One  can  see,  in  particular,  a  decrease  of 


Fig.  12.  Distributions  of  the  current  densities  along  the  SOFC(Q2  ,  H+)  for  the  co-flow  (thin  lines)  and  counter-flow  (bold  lines)  modes.  Case  of  pr  =  0.7. 


Table  1 

The  main  characteristics  of  the  SOFC(02-,  H+)  in  comparison  with  the  SOFC(02~)  and  the  SOFC(H+) 


7h 

Feeding  mode 

Case  of  the  maximum  voltage 

Case  of  pr  =  0.7 

Peff  P av  (mV) 

m 

U  (mV) 

1 

Peff  P av  (mV) 

m 

U  (mV) 

r\ 

0 

Co-flow 

88.1 

0.9 

908 

0.650 

333.6 

0.900 

111 

0.556 

0.5 

Co-flow 

91.0 

%1 

960 

0.766 

393.2 

0.954 

821 

0.625 

Counter-flow 

94.0 

%1 

980 

0.782 

394.4 

0.954 

822 

0.626 

1 

Co-flow 

201.8 

%1 

1012 

0.807 

434.5 

0.998 

865 

0.689 

Counter-flow 

123.2 

%1 

1057 

0.842 

436.0 

0.998 

865 

0.689 
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the  terminal  voltage  leads  to  a  decrease  of  the  fuel  utilization 
(Fig.  10),  to  an  increase  of  the  emf  values  (Fig.  11),  and  to 
an  increase  of  the  current  density  values  and  to  their  higher 
uniformity  (Fig.  12). 

The  most  important  characteristics  of  the  SOFC,  such  as 
the  terminal  voltage,  the  average  specific  power  calculated 
according  Eq.  (14),  the  fuel  utilization,  and  the  efficiency 
are  presented  in  the  Table  1  for  the  S0FC(02~,  H+)  in 
comparison  with  the  S0FC(02-)  and  the  SOFC(H+).  One 
can  see  that  the  higher  proton  transfer  number  the  higher  the 
SOFC  efficiency.  Increase  of  the  terminal  voltage  also  leads 
to  a  noticeable  increase  of  the  SOFC  efficiency,  however, 
this  leads  to  a  considerable  decrease  of  the  SOFC  specific 
power. 

4.  Conclusions 

The  analysis  presented  shows  that  the  SOFC  based  on 
the  co-ionic  electrolyte  has  significant  advantages  over 
the  SOFC  based  on  the  oxygen  ion  electrolyte  and  very 


close  to  the  efficiency  of  the  SOFC  based  on  the  pro¬ 
ton  electrolyte  when  they  operate  at  the  relative  power  of 
0.7.  The  S0FC(02-,  H+)  efficiency  is  very  high  when 
it  operates  at  maximum  power,  but  the  specific  power  in 
this  case  is  significantly  lower.  The  main  aim  in  SOFC 
technology  must  be  the  elaboration  of  a  cheap  and  stable 
co-ionic  electrolyte  with  as  high  proton  transfer  number  as 
possible. 
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